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In this paper, we present results from a simulation of stellar core collapse, bounce, and postbounce evolution with 
Boltzmann neutrino transport. We motivate the development of our Boltzmann solver in light of the sensitivity 
of the neutrino-heating core collapse supernova paradigm to details in the neutrino transport, particularly near 
the neutrinospheres, where the neutrinos are neither diffusing nor free streaming and a kinetic description is 
necessary, and in light of the mixed outcomes and transport approximations used in all prior supernova models 
in both one and two dimensions. We discuss the implications of our findings for the supernova mechanism and 
future supernova research. We also present the results of a Boltzmann transport prediction of the early neutrino 
light curves in the model included here. 



1. Introduction 

Beginning with the first numerical simulations con- 
ducted by Colgate and White jjj], three decades of 
supernova modeling have established a basic super- 
nova paradigm. The supernova shock wave — formed 
when the iron core of a massive star collapses grav- 
itationally and rebounds as the core matter exceeds 
nuclear densities and becomes incompressible — stalls 
in the iron core as a result of enervating losses to nu- 
clear dissociation and neutrinos. The failure of this 
"prompt" supernova mechanism sets the stage for a 
"delayed" mechanism, whereby the shock is reener- 
gized by the intense neutrino flux emerging from the 
neutrinospheres carrying off the binding energy of the 
proto- neutron starj^j^]. The heating is mediated pri- 
marily by the absorption of electron neutrinos and 
antineutrinos on the dissociation-liberated nucleons 
behind the shock. This past decade has also seen 
the emergence of multidimensional supernova models, 
which have investigated the role convection, rotation, 
and magnetic fields may play in the explosion [p|-|lot, 
in some cases invoking new explosion paradigms. 

Although a plausible framework is now in place, 
fundamental questions about the explosion mecha- 
nism remain: Is the neutrino heating sufficient, or 
are multidimensional effects such as convection and 



rotation necessary? Can the basic supernova observ- 
able, explosion, be reproduced by detailed spherically 
symmetric models, or are multidimensional models 
required? Without a doubt, core collapse supernovae 
are not spherically symmetric. For example, neu- 
tron star kicks [ h l | and the polarization of supernova 
emitted light |1 2| cannot arise in spherical symmetry. 
Nonetheless, ascertaining the explosion mechanism 
and understanding every explosion observable are two 
different goals. To achieve both, simulations in one, 
two, and three dimensions must be coordinated. 

2. Convection 

Supernova convection falls into two categories: (1) 
convection near or below the neutrinospheres, which 
we refer to as proto-neutron star convection and (2) 
convection between the gain radius and the shock, 
which we refer to as neutrino-driven convection. The 
gain radius is the radius at which neutrino heating 
and cooling via electron neutrino and antineutrino 
absorption and emission between the neutrinospheres 
and the shock balance. There is net neutrino heating 
above this radius and net neutrino cooling below it. 

Proto-neutron star convection may aid the explo- 
sion mechanism by boosting the neutrinosphere lumi- 
nosities. Hot, lepton-rich rich matter is convectively 
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transported to the neutrinospheres. This mode of 
convection may develop owing to instabilities caused 
by lepton and entropy gradients established by the 
deleptonization of the proto-neutron star via elec- 
tron neutrino escape near the electron neutrinosphere 
and by the weakening supernova shock (as the shock 
weakens, it causes a smaller entropy jump in the ma- 
terial flowing through it). Proto-neutron star convec- 
tion is arguably the most difficult to investigate nu- 
merically because the neutrinos and the matter are 
coupled and, consequently, multidimensional simula- 
tions must include both multidimensional hydrody- 
namics and multidimensional, multigroup neutrino 
transport. [Multigroup, i.e., multi-neutrino energy, 
transport is necessary because the neutrino opacities 
are strongly energy dependent and low- and high- 
energy neutrinos may be transported in very differ- 
ent ways (e.g., diffusion versus free streaming) at any 
given spatial point in the core at any given time.] 

Neutrino-driven convection may aid the explosion 
mechanism by boosting the shock radius and the neu- 
trino heating efficiency, thereby facilitating shock re- 
vival, ft develops as the result of the entropy gradient 
established as the shocked stellar core material infalls 
between the shock and the gain radius, being contin- 
ually heated in the process. 

2.1. Proto-Neutron Star Convection (ID): 
Neutron Fingers 

The fundamental difficulty in modeling convection 
in spherically symmetric models is apparent: convec- 
tion is a three-dimensional phenomenon, and spher- 
ically symmetric models can incorporate convection 
only in a phenomenological way (e.g., via a mixing- 
length prescription). Moreover, because convection is 
not admitted by the one-dimensional hydrodynamics 
equations, some imposed criterion for the existence of 
convection must be used. 

Neutron-finger convection has been invoked by 
Wilson et al.p3| in their one-dimensional models and 
has been deemed necessary by them to obtain su- 
pernova explosions. This mode of proto-neutron star 
convection arises in the presence of a negative electron 
fraction gradient and a positive entropy gradient in 
the postshock stellar core, resulting in higher-entropy, 
neutron-richer matter above lower-entropy, neutron- 
poorer matter in the core. With the assumption that 
energy transport by neutrinos is more efficient than 
lepton transport, neutron fingers develop under these 
conditions, resulting (like salt fingers in the ocean) in 
finger-like downflows of neutron-rich matter that pen- 
etrate deep into the stellar core. The assumption that 



energy transport is more efficient than lepton trans- 
port is justified in the following way: Three flavors 
of neutrinos (electron, muon, and tau) can transport 
energy, whereas only one (electron) can transport lep- 
ton number. However, detailed neutrino equilibration 
experiments carried out by Bruenn and Dineva|l4j 
demonstrate that the muon and tau neutrinos do not 
couple strongly with the stellar core matter in energy, 
and therefore, there is only one flavor (electron) that 
transports both energy and lepton number efficiently. 
Given the outcome of these numerical experiments, 
the fundamental assumption made by Wilson et al. 
should be reexamined. However, it is also important 
to note that the equilibration experiments carried out 
by Bruenn and Dineva have to be repeated in light 
of energy exchange channels between the muon and 
tau neutrinos and the stellar core matter that have 
recently been identified[fl5[. 

2.2. Proto-Neutron Star Convection (2D): 
Ledoux Convection 

In certain regions of the stellar core, neutrino trans- 
port can equilibrate a convecting fluid element with 
its surroundings in both entropy and lepton number 
on time scales shorter than convection time scales, 
rendering the fluid element nonbouyant. This will oc- 
cur in intermediate regimes in which neutrino trans- 
port is efficient but in which the neutrinos are still 
strongly enough coupled to the matter. Figures 1 
and 2 from Mezzacappa et al.[j?j] demonstrate that 
this equilibration can in fact occur. Figure 1 shows 
the onset and development of proto-neutron star con- 
vection in a 25 Mq model shortly after bounce in a 
simulation that did not include neutrino transport, 
i.e., that was a hydrodynamics-only run. Figure 2 on 
the other hand shows the lack of any significant onset 
and development of convection when neutrino trans- 
port was included in what was otherwise an identi- 
cal model. Transport's damping effects are obvious. 
(The same result occurred in our 15 M model 71) 

On the other hand, in the model of Keil et ah] lq] , 
vigorous proto-neutron star convection developed, 
which then extended deep into the core as a delep- 
tonization wave moved inward owing to neutrinos dif- 
fusing outward. In this model, convection occurs very 
deep in the core where neutrino opacities are high and 
transport becomes inefficient in equilibrating a fluid 
element with its surroundings. 

It is also important to note in this context that 
Mezzacappa et al. and Keil et al. used comple- 
mentary transport approximations. In the former 
case, spherically symmetric transport was used, which 
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Figure 1. Two-dimensional entropy plots showing 
the evolution of proto-neutron star convection in our 
hydrodynamics-only 25 Mq model at 12, 17, and 27 
ms after bounce. 
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2.3. Neutrino-Driven Convection (2D) 

This mode of convection occurs directly between 
the gain radius and the stalled shock as a result of 
the entropy gradient that forms as material infalls 
between the two while being continually heated from 
below. In Figure 3, a sequence of two-dimensional 
plots of entropy are shown, illustrating the develop- 
ment and evolution of neutrino-driven convection in 
our 15 Mq model Ji|. High-entropy, rising plumes and 
lower-entropy, denser, finger-like downflows are seen. 
The shock is distorted by this convective activity. 




Figure 3. Two-dimensional entropy plots showing 
the evolution of neutrino-driven convection in our 15 
Mq model at 137, 212, and 512 ms after bounce. 



Figure 2. Two-dimensional entropy plots showing 
the evolution of proto-neutron star convection in our 
hydrodynamics-plus-neutrino-transport 25 Mq model 
at 12, 17, and 27 ms after bounce. 



maximizes lateral neutrino transport and overesti- 
mates the neutrino-matter equilibration rate; in the 
latter case, ray-by-ray transport was used, which min- 
imizes (zeroes) lateral transport and underestimates 
the neutrino-matter equilibration rate. 

These outcomes clearly demonstrate that to deter- 
mine whether or not proto-neutron star convection 
exists and, if it exists, is vigorous, will require sim- 
ulations coupling three-dimensional, multigroup neu- 
trino transport and three-dimensional hydrodynam- 
ics. Moreover, realistic high-density neutrino opaci- 
ties will also be needed. 



In the Herant et al.|J simulations, large-scale con- 
vection developed beneath the shock, leading to in- 
creased neutrino energy deposition, the accumulation 
of mass and energy in the gain region, and a ther- 
modynamic engine they claimed ensured explosion, 
although Herant et al. stressed the need for more 
sophisticated multidimensional, multigroup transport 
in future models. [They used two-dimensional "gray" 
(neutrino-energy-integrated, as opposed to multi- 
group) flux-limited diffusion in neutrino-thick regions 
and a neutrino lightbulb approximation in neutrino- 
thin regions. In a lightbulb approximation, the neu- 
trino luminosities and rms energies are assumed con- 
stant with radius.] In the Burrows et al. simu- 
lations ||, neutrino-driven convection in some mod- 
els significantly boosted the shock radius and led to 
explosions. However, they stressed that success or 
failure in producing explosions was ultimately deter- 
mined by the values chosen for the neutrino spectral 
parameters in their gray ray-by-ray (one-dimensional) 
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neutrino diffusion scheme. [In spherical symmetry 
(ID), all rays are the same. In a ray-by-ray scheme 
in axisymmetry (2D), not all rays are the same, al- 
though the transport along each ray is a ID prob- 
lem. In this latter case, lateral transport between 
rays is ignored.] Focusing on the neutrino lumi- 
nosities, Janka and MullerQ, using an adjustable 
central neutrino lightbulb, conducted a parameter 
survey and concluded that neutrino-driven convec- 
tion aids explosion only in a narrow luminosity win- 
dow (±10%), below which the luminosities are too 
low to power explosions and above which neutrino- 
driven convection is not necessary to power explo- 
sions. In more recent simulations carried out by 
Swestyjnj] using two-dimensional gray flux-limited 
diffusion in both neutrino-thick and neutrino-thin re- 
gions, it was demonstrated that the simulation out- 
come varied dramatically as the matter-neutrino "de- 
coupling point," which in turn sets the neutrino spec- 
tra in the heating region, was varied within reason- 
able limits. (The fundamental problem in gray trans- 
port schemes is that the neutrino spectra, which are 
needed for the heating rate, are not computed. The 
spectra are specified by choosing a neutrino "tem- 
perature," normally chosen to be the matter tem- 
perature at decoupling. In a multigroup scheme, 
the spectra are computed dynamically.) In our two- 
dimensional models M, the angle-averaged shock radii 
do not differ significantly from the shock trajecto- 
ries in their one-dimensional counterparts, and no 
explosions are obtained, as seen in Figure 3. Nei- 
ther the luminosities nor the neutrino spectra are free 
parameters. Our two-dimensional simulations im- 
plemented precomputed spherically symmetric (ID) 
multigroup flux-limited diffusion neutrino transport, 
compromising transport dimensionality to implement 
multigroup transport and a seamless transition be- 
tween neutrino-thick and neutrino-thin regions, al- 
though without feedback between the hydrodynamics 
and the transport. 

In light of the neutrino transport approximations 
made, the fact that none of the simulations have 
been three dimensional, and the mixed outcomes, 
next-generation simulations will have to reexplore 
neutrino-driven convection in the context of three- 
dimensional simulations that implement more realis- 
tic three-dimensional multigroup neutrino transport. 



3. General Relativity, Rotation, 
and Magnetic Fields 

For discussions of the role of general relativity, 
rotation, and magnetic fields in supernova models, 
the reader may begin with the papers by Bruenn 
et al.lOJ, Liebendorfer et al.[Q,|20), Fryer and 
HegerQ7 Khokhlov et al.|loj, and MacFadyen and 
WoosleyQ. 



4. Boltzmann Neutrino Transport (ID) 

The neutrino energy deposition behind the shock 
depends sensitively not only on the neutrino lumi- 
nosities but also on the neutrino spectra and angular 
distributions in the postshock region, necessitating 
exact multigroup Boltzmann neutrino transport near 
and above the neutrinospheres or a very good approx- 
imation of it. Ten percent variations in any of these 
quantities can make the difference between explosion 
and failure in supernova models Past simu- 

lations have implemented increasingly sophisticated 
approximations to multigroup Boltzmann transport, 
the most sophisticated of which is multigroup flux- 
limited diffusion A generic feature of this ap- 
proximation is that it underestimates the isotropy of 
the neutrino angular distributions in the heating re- 
gion and, thus, the heating rate[ p4p5| ]. Therefore, 
the question arises whether or not failures to produce 
explosions in past one-dimensional models were the 
result of the transport approximations employed. It 
is important to note that, without invoking proto- 
neutron star (e.g., neutron finger) convection, sim- 
ulations that implement multigroup flux-limited dif- 
fusion do not produce explosions p3|Jl3] ] (as we will 
discuss, the existence and vigor of proto-neutron star 
convection is a matter of debate fj], 14|Jl^ ] ) . 

To begin to address the question posed above, we 
have been simulating the core collapse, bounce, and 
postbounce evolution of 13, 15, and 20 Mq stars, 
beginning with the precollapse models of Nomoto 
and Hashimoto^fjj, with a new neutrino radiation 
hydrodynamics code for both Newtonian and gen- 
eral relativistic spherically symmetric flows: AGILE- 
BOLTZTRAN. BOLTZTRAN is a three-flavor Boltz- 
mann neutrino transport solver [E7 U 2 8 ] , now extended 
to fully general relativistic flowsfll9[]T In the 13 Mq 
Newtonian (gravity) simulation we present here, it is 
employed in the 0(v/c) limit. AGILE is a conserva- 
tive, adaptive mesh, general relativistic hydrodynam- 
ics code]l^,^|. Its adaptivity enables us to resolve 
and seamlessly follow the shock through the iron core 
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into the outer stellar layers. 

The equation of state of Lattimer and Swesty|3o| 
(LS EOS) is employed to calculate the local thermo- 
dynamic state and nuclear composition of the mat- 
ter in nuclear statistical equilibrium (NSE). For mat- 
ter initially in the silicon layer, the temperatures 
are insufficient to achieve NSE. In this region, the 
radiation and electron components of the LS EOS 
are used, while an ideal gas of 28 Si is assumed for 
the nuclear component. For typical hydrodynamic 
timesteps (~ .1 millisecond), silicon burning occurs 
within a single timestep for T ~ 5GK[^I); there- 
fore, when a fluid element exceeds a temperature of 5 
GK in our simulation, the silicon is instantaneously 
burned, achieving NSE and releasing thermal energy 
equal to the difference in nuclear binding energy be- 
tween 28 Si and the composition determined by the LS 
EOS. 

Figure 4, taken from the simulation of Mezzacappa 
et al.|3^], shows the radius- versus-time trajectories of 
equal mass (O.OIM0) shells in the stellar iron core and 
silicon layer in a Newtonian simulation initiated from 
the 13 Mq progenitor. Core bounce and the forma- 
tion and propagation of the initial bounce shock are 
evident. This shock becomes an accretion shock, de- 
celerating the core material passing through it. At 
ps 125 ms after bounce, the accretion shock stalls at 
a radius w 250 km and begins to recede, continuing 
to do so during the first 500 ms of postbounce evolu- 
tion. No explosion has developed in this model dur- 
ing this time. Similar behavior is exhibited in our 15 
and 20 Mq Newtonian models and in our 13 and 20 
Mq general relativistic models, although these have 
not yet reached 500 ms after bounce. We will con- 
tinue these simulations and report on the final out- 
comes in subsequent papers. If explosions are consis- 
tently obtained in the more realistic general relativis- 
tic cases, this would imply that core collapse super- 
novae driven purely by neutrino heating are possible, 
although these same models would have to be con- 
sidered in three dimensions to assess whether or not 
multidimensional effects alter this conclusion and, of 
course, to model any associated phenomenology, such 
as neutron star kicks. 

Figure 5 shows the time evolution of the three- 
flavor neutrino signal computed with Boltzmann neu- 
trino transport shortly after shock breakout in our 
general relativistic simulation^]. We see the elec- 
tron neutrino burst and the three-flavor emission de- 
velop from the hot, shocked mantle. This early evolu- 
tion is a consequence of the time-dependent neutrino 
transport in semitransparent regions. 
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Figure 4. We trace the shock, nuclear burning, and 
dissociation fronts (the shock and dissociation fronts 
are coincident), which carve out three regions in the 
(r, t) plane. A: Silicon. B: Iron produced by infall 
compression and heating. C: Free nucleons and alpha 
particles. 



5. Outlook 

Fundamental questions in supernova theory re- 
main, which can only be answered via systematic and 
coordinated simulations in one, two, and three dimen- 
sions. 

We have shown results from the first 500 ms of our 
one-dimensional (spherically symmetric) Newtonian 
simulation with Boltzmann neutrino transport initi- 
ated from a 13 Mq progenitor. In light of our imple- 
mentation of Boltzmann transport, if we do not ob- 
tain explosions in this model, or other models we have 
initiated from different progenitors (see also Rampp 
and Janka|Q), it would suggest that improvements 
in our initial conditions (precollapse models) and/or 
input physics are needed, and/or that the inclusion 
of multidimensional effects such as convection, rota- 
tion, and magnetic fields are required ingredients in 
the recipe for explosion. In the past, it was not clear 
whether failure in spherically symmetric models was 
the result of transport approximations or the lack of 
inclusion of important physics. With the implemen- 
tation of Boltzmann transport, this conclusion can be 
made unambiguously. We will report on the contin- 
ued evolution of our 13 Mq model and on our other 
models in subsequent papers. 
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convection in the proto-neutron star, which in turn 
will depend on the neutrino transport, all of these 
studies must be strongly coupled. 

We have developed a general relativistic neutrino 
Boltzmann transport /radiation hydrodynamics code, 
AGILE-BOLTZTRAN, that can be used to study the 
supernova mechanism and nucleosynthesis, and to 
make accurate predictions of the neutrino signatures 
in supernovae and failed supernovae. In a model initi- 
ated from a 13 Mq progenitor, we have computed the 
early three-flavor neutrino signal, with general rela- 
tivistic Boltzmann neutrino transport. We are cur- 
rently running other models with different progenitor 
masses and will report on their dynamics and neu- 
trino signatures in future papers. 



Figure 5. We plot the three- flavor neutrino luminosi- 
ties and rms energies at 500 km as a function of time 
over the first 40 ms after bounce. 



Potential improvements in our initial conditions 
and input physics include: improvements in precol- 
lapse models [M-pTj; the use of ensembles of nuclei in 
the stellar core rather than a single representative nu- 
cleus; computing the neutrino-nucleus cross sections 
with detailed shell model computations and the 
inclusion of nucleon correlations in the high-density 
neutrino opacities |i9|JIo] ] . These improvements all 
have the potential to quantitatively, if not qualita- 
tively, change the details of our simulations. Thus, it 
is important to note that the conclusions drawn here 
are drawn considering the initial conditions and input 
physics used. 

To accurately investigate multdimensional effects 
such as convection, rotation, and magnetic fields, 
future simulations must be carried out in three 
dimensions and must implement realistic, three- 
dimensional, multigroup neutrino transport. Three- 
dimensional simulations will be necessary to assess, 
for example, the vigor of convection in the proto- 
neutron star, where the neutrinos and the matter are 
strongly coupled and the flow is three-dimensional, 
and to assess the character of neutrino-driven convec- 
tion behind the shock in a stellar core that is both ro- 
tating and convecting. Certainly, three-dimensional 
simulations are required to study the development of 
MHD jets in stellar cores, and given that the devel- 
opment of such jets depends in some scenarios on the 
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